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Abstract. Results from a model study on the formation of exotic quark clusters 
at the hadronization stage of a heavy ion collision are presented. The dynamical 
quark molecular dynamics (qMD) model which is used is sketched, and results 
for exotica made of up to six (anti-)quarks are shown. The second part focuses 
on pentaquarks. The rapidity distribution are shown, and the distribution of 
strangeness is found to yield an indicator of thermalization and homogcnisation 
of the deconfined quark system. Relative 0+ yields are found to be lower than 
thermal model estimates. 
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In the quark model of QCD, quarks are associated with fundamental representations 
of the colour-gauge group SU(3). The low energy regime of QCD contains only colour 
singlets such as mesons and baryons, which are made up from a quark and anti-quark 
with colour and anti-colour, or three quarks with three different colours, respectively. 
However, colour neutrality can be fulfilled whenever the net number of quarks is a 
multiple of three, thus opening the possibility of so called exotic hadrons such as the 
tetraquark (or dimeson), the pentaquark, and so on. New interest in the study of 
exotic hadrons has been created by recent reports about the discovery of pentaquark 
states 8 + and 5 in pA, nA, and 7A collisions (P] - see for discussions of the 
theoretical understanding) . 

Here, we will discuss the formation of clusters of up to six quarks or anti-quarks 
in the hadronization stage of a heavy ion collision, with a focus on pentaquarks. In 
such an event, two nuclei collide, creating hot and dense matter, which can set free 
colour charges if the energy density is high enough. During the expansion of the hot 
system, colour charges will recombine to form colour neutral hadrons which finally 
can be detected. The study is motivated by the idea that the hadronization from a 
hot quark soup is an environment favourable to the formation of exotica [Hj • 

In order to describe the dynamics of the collision, the hadronic model UrQMD 
is used to create an initial state of hot and dense matter. The particle content 
after complete overlap of the nuclei is converted to quarks, and the quark molecular 
dynamics model qMD is used to describe the subsequent dynamics and recombination 
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Figure 1. Rapidity distribution of different types of quark clusters at upper 
SPS energies (left figure). Pentaquark clusters are represented by the thick doted 
line. Rapidity distribution of pentaquark clusters, normalised to the total baryon 
distribution, for different colliding systems (right figure). 



of the quark system 4 . In the qMD model, quarks are described as classical particles 
with current masses subject to the Hamiltonian 

i— 1 i,j \ \ 1 J \ / 

All quarks carry spin and isospin which is used in the subsequent mapping to hadrons, 
and colour charges Wi represented by three vectors in the root diagram, with anti- 
quarks carrying the corresponding anti-colours. No colour exchange between quarks 
is assumed, leaving the colour of a particle unchanged throughout the complete time 
evolution of the system. The sign of the potential interaction between two particles 
depends on the scalar product of the root vectors, 

CP- = wjw 3 , (2) 

yielding repulsion for equal charges and attraction between charge and anti-charge. 
The shape of the potential interaction is motivated by the Cornell potential, where in 
all calculations the Coulomb term is omitted, and the string constant k is treated as a 
parameter of the model. The time evolution yields an expanding system where quarks 
form colour neutral clusters. Once a cluster is separated in space from the remaining 
quarks and its interaction with the system falls below a threshold, it is mapped to a 
hadron, taking into account the quantum numbers and energy-momentum four-vectors 
of the quarks in the cluster |S] . 

We investigate the formation of colour neutral clusters with quark configurations 
of mesons (qq), baryons (qqq), dimesons (qqqq), pentaquarks (qqqqq), dibaryons 
(qqqqqq), baryonia {qqqqqq), and the corresponding anti-particles. We cannot make 
any statement about the actual existence and physical properties of exotic hadrons 
from our extremely simplistic classical model, but we assume that any physical exotic 
hadron formed by quark coalescence at hadronization should emerge from a colour 
neutral cluster as described in the qMD model. The dynamical evolution of systems 
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Figure 2. Distribution over strangeness of pentaquark cluster yield (a), and 
normalised to the combinatorial expectation value (b). Yields of © + (S = +1 with 
a s quark) and 3 (S = —2 with a q quark) are only half of the combinatorial 
value at lower SPS energies. 



at upper SPS energies then yields the rapidity distribution of clusters as shown on the 
left part of figure ^ 

Focusing on pentaquark clusters, the right part of figure ^ shows the rapidity 
distribution normalised to the distribution of baryons, which is relevant for the chances 
of identifying pentaquarks in the debris of a heavy ion collision. At least for SPS 
energies, this ratio is nearly twice as large at Ay = ±1 off midrapidity as at y = 0. 
For RHIC energies, a cylindrical quark gas, thermalised at a temperature of 250 MeV, 
was assumed as initial condition, so that the full dynamics of the expanding fireball 
is missing. The apparently high relative yield of pentaquarks is strongly reduced, of 
course, by the large size of the flavour multiplet, 

3®3®3(g>3®3 = 3- l©8-8®4-10©2-l0©3-27e35, (3) 

where, assuming SU(3) coupling, only 1 out of 243 states corresponds to the 9 + , while 
the clustering procedure populates the whole multiplet structure. The left part of 
figure [21 shows the relative distribution of pentaquark clusters over strangeness. The 
grey band corresponds to the fully SU(3) symmetric case. For SPS energies, there 
are large deviations, corresponding to the strangeness suppression. But even when 
normalising the relative distribution to the actual strangeness content of the quark 
system (right part of figureEJ), deviations remain for SPS energies. In our model, this 
is due to the correlations carried over into the quark system from the original hadrons. 
However, we conclude that the relative strangeness of pentaquarks can give a measure 
for thermalization and homogenisation of the deconfined quark system. 

Finally, we compare the relative yield of + / A as a function of bombarding energy 
from our quark clustering dynamics to the results obtained from grand canonical 
estimates [H] and thermal models including strangeness suppression factors [7] 
(Figure El. We find relative yields roughly one order of magnitude below the earlier 
estimates. The large downward range is due to the use of different coupling schemes, 
SU(2) for fixed strangeness, SU(3), and SU(6) for spin and flavour. It should be noted 



Exotic hadrons from dynamical clustering of quarks 



4 



g 
CO 



1 



10' 



o 

♦ 



gee (Randrup) 

thermal fit (Letessier et £ 

qMD 



10' 



10" 



/s [GeV] 



Figure 3. Ratio 0+/A as a function of energy from the grand canonical 
ensemble, thermal models including strangeness suppression factors, and from 
qMD clustering. 



that the qMD result does not include hadronic rescattering, which will both destroy 
<d + s formed from quark clustering and create new ones from hadronic coalescence 
(HU - there are, however, much unknowns about the formation of pentaquarks, see e. 
g- PU), nor the possible population of + from decaying N* pentaquarks (see e. g. 0). 
These factors shift the + /A and require further study. 
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